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The quantum ferromagnetic transition of itinerant electrons is considered. We give a pedagogical 
review of recent results which show that zero-temperature soft modes that are commonly neglected, 
invalidate the standard Landau- Ginzburg- Wilson description of this transition. If these modes are 
taken into account, then the resulting order parameter field theory is nonlocal in space and time. 
Nevertheless, for both disordered and clean systems the critical behavior has been exactly deter- 
mined for spatial dimensions d > 2 and d > 1, respectively. The critical exponents characterizing the 
paramagnetic-to-ferromagnetic transition are dimensionality dependent, and substantially different 
from both mean-field critical exponents, and from the classical Heisenberg exponents that charac- 
terize the transition at finite temperatures. Our results should be easily observable, particularly 
those for the disordered case, and experiments to check our predictions are proposed. 



I. INTRODUCTION 

Phase transitions that occur in a quantum mechani- 
cal system at zero temperature (T = 0) as a function of 
some non-thermal control parameter are called quantum 
phase transitions. In contrast to their finite-temperature 
counterparts, which are often referred to as thermal or 
classical phase transitions, the critical fluctuations one 
has to deal with at zero temperature are quantum fluc- 
tuations rather than thermal ones, and the need for a 
quantum mechanical treatment of the relevant statistical 
mechanics makes the theoretical description of quantum 
phase transitions somewhat different from that of clas- 
sical ones. However, as Hertz has shown in a seminal 
paperjj the basic theoretical concepts that have led to 
successfully describe and understand thermal transitions 
work in the quantum case as well. 

Experimentally, the zero-temperature behavior of any 
material can of course not be studied directly, and fur- 
thermore the most obvious control parameter that drives 
a system through a quantum transition is often some mi- 
croscopic coupling strength that is hard to change exper- 
imentally. As a result, the dimensionless distance from 
the critical point, t, which for classical transitions with a 
transition temperature T c is given by t = T/T c — 1 and is 
easy to tune with high accuracy, is much harder to con- 
trol in the quantum case. However, t is usually dependent 
on some quantity that can be experimentally controlled, 
like e.g. the composition of the material. Also, the zero 
temperature critical behavior manifests itself already at 
low but finite temperatures. Indeed, in a system with a 
very low thermal transition temperature all but the final 
asymptotic behavior in the critical region is dominated 



by quantum effects. The study of quantum phase tran- 
sitions is therefore far from being of theoretical interest 
only. 

Perhaps the most obvious example of a quantum phase 
transition is the paramagnet-to-ferromagnet transition of 
itinerant electrons at T = as a function of the exchange 
interaction between the electronic spins. Early theoret- 
ical workEl on this transition suggested that the critical 
behavior in the physical dimensions d = 2 and d = 3 was 
not dominated by fluctuations, and mean-field like, as is 
the thermal ferromagnetic transition in dimensions d > 4. 
The reason for this is a fundamental feature of quantum 
statistical mechanics, namely the fact that statics and 
dynamics are coupled. As a result, a quantum mechan- 
ical system in d dimensions is very similar so the corre- 
sponding classical system in d + z dimensions, where the 
so-called dynamical critical exponent z can be thought of 
as an extra dimensionality that is provided to the system 
by time or temperature. The d + z-dimensional space rel- 
evant for the statistical mechanics of the quantum system 
bears some resemblance to d + 1-dimensional Minkowski 
space, but z does not need to be equal to 1 in nonrelativis- 
tic systems. For clean and disordered itinerant quantum 
ferromagnets, one finds z — 3 and z = 4, respectively, 
in mean-field theory. This appears to reduce the up- 
per critical dimension above which fluctuations are 
unimportant and simple mean-field theory yields the cor- 
rect critical behavior, from rf+ = 4 in the classical case 
to dj = 1 and d+ = 0, respectively, in the clean and 
disordered quantum cases. If this were true, then this 
quantum phase transition would be rather uninteresting 
from a critical phenomena point of view. 

It has been known for some time that, for the case of 



1 



disordered systems, this conclusion cannot be correct. El 
It is known that in any system with quenched disor- 
der that undergoes a phase transition, the critical ex- 
ponent v that describes the divergence of the correlation 
length, t~ v for t — ► 0, must satisfy the inequality 
v > 2/d.cl However, mean-field theory yields v = 1/2, 
which is incompatible with this inequality for d < 4. 
Technically, this implies that the disorder must be a rel- 
evant perturbation with respect to the mean-field fixed 
point. The mean- field fixed point must therefore be un- 
stable, and the phase transition must be governed by 
some other fixed point that has a correlation length ex- 
ponent v > 2/d. 

Recently such a non-mean field like fixed point has 
been discovered, and the critical behaviot,has been de- 
termined exactly for all dimensions d > 2.Q It was found 
that both the value c?+ = for the upper critical dimen- 
sion, and the prediction of mean-field critical behavior for 
d > <i+ were incorrect. Instead, <f+ = 2, and while both 
the quantum fluctuations and the disorder fluctuations 
are irrelevant with respect to the new fixed point for all 
d > e?+ , there are two other "upper critical dimensional- 
ities", d^ + = 4 and d+ ++ = 6. The critical behavior for 
e?+ < d < <i+ ++ is governed by a non-standard Gaussian 
fixed point with non-mean field like exponents, and only 
for d > d^ ++ does one obtain mean-field exponents. In 
addition, the clarification of the physics behind this sur- 
prising behavior has led to the conclusion that very sim- 
ilar effects occur in clean systems.Q In that case, d+ = 1 
in agreement with the early result, but again the criti- 
cal behavior is nontrivial in a range of dimensionalities 
d+ < d < d++ = 3, and only for d > d++ does one obtain 
mean-field critical behavior. In addition, we have found 
that Hertz's 1 — e expansion for the clean case is invalid. 
This explains an inconsistency between this expansion 
and an exact exponent relation that was noted earlier 
by Sachdev.Q In order to keep our discussion focused, in 
what follows we will restrict ourselves to the disordered 
case, where the effects are more pronounced, and will 
only quote results for the clean case where appropriate. 

The basic physical reason behind the complicated be- 
havior above the upper critical dimensionality <i+, i.e. in 
a regime in parameter space where the critical behavior 
is not dominated by fluctuations, is simple. According 
to our general understanding of continuous phase transi- 
tions or critical points, in order to understand the critical 
singularities at any such transition, one must identify all 
of the slow or soft modes near the critical point, and one 
must make sure that all of these soft modes are properly 
included in the effective theory for the phase transition. 
This is obvious, since critical phenomena arc effects that 
occur on very large length and time scales, and hence 
soft modes, whose excitation energies vanish in the limit 
of long wavelengths and small frequencies, will in gen- 
eral influence the critical behavior. In previous work on 
the ferromagnetic transition it was implicitly assumed 
that the only relevant soft modes are the fluctuations 
of the order parameter, i.e. the magnetization. For fi- 



nite temperatures this is correct. However, at T = 
there are additional soft modes in a disordered electron 
system, namely diffusive particle-hole excitations that 
arc distinct from the spin density excitations that de- 
termine the magnetization. In many-body perturbation 
theory these modes manifest themselves as products of 
retarded and advanced Green's functions, and in field 
theory they can be interpreted as the Goldstone modes 
that result from the spontaneous breaking of the symme- 
try between retarded and advanced correlation functions, 
or between positive and negative imaginary frequencies. 
In a different context, namely the transport theory for 
disordered electron systems, these diffusive excitations 
are sometimes referred to as 'diffusons' and 'Cooper- 
ons', respectively, and they are responsible for what is 
known as 'weak localization effects' in disordered elec- 
tron systems.0 For our purposes, their most important 
feature is their spatial long-range nature in the zero fre- 
quency limit. This long-range nature follows immediately 
from the diffusion equation 

{d t -Ddl) /(x,f) = , (1.1a) 

for some diffusive quantity /, with D the diffusion con- 
stant. Solving this equation by means of a Fourier- 
Laplace transform to wavevectors q and complex frequen- 
cies z, one obtains in the limit of zero frequency, 

/(q,« = 0) = -L/(q,t = 0) . (1.1b) 

Long-range static correlations are thus an immediate con- 
sequence of the diffusive nature of the density dynamics 
in disordered systems. 

The fact that we are concerned with the zero frequency 
or long-time limit is due to the order parameter, i.e. the 
magnetization, being a conserved quantity. Since the 
only way to locally change the order parameter density 
is to transport this conserved quantity from one region 
in space to another, in order to develop long-range order 
over arbitrarily large distances the systems needs an in- 
finitely long time. This in turn means that criticality can 
be reached only if the frequency is taken to zero before 
the wavenumbcr. This feature would be lost if there were 
some type of spin-flip scattering mechanism present, and 
our results hold only in the absence of such processes. 
For the same reason, they do not apply to quantum an- 
tiferromagnets, which show a quite different behavioral 

It is important that the long-range static correlations 
mentioned above are distinct from the order parameter 
fluctuations. For instance, the latter are soft only at the 
critical point and in the ordered phase, while the former 
are soft even in the paramagnetic phase, and they do 
not change their nature at the critical point. However, 
since they couple to the conserved order paramter, they 
influence the critical behavior. If one integrates out these 
diffusive modes in order to obtain an effective theory or 
Landau-Ginzburg- Wilson (LGW) functional in terms of 
the order parameter only, then their long-range nature 
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leads to infrared singular integrals, which in turn results 
in singular vertices in the LGW funcional, or diverging 
coupling constants for couplings between the order pa- 
rameter fluctuations. The usual LGW philosophy of de- 
riving an effective local field theory entirely in terms of 
the order parameter field therefore does not lead to a well 
behaved field theory in this case. The situation is analo- 
gous to a well known phenomenon in high energy physics: 
Suppose some interaction between, say, fermions, is me- 
diated by the exchange of some other particles, e.g. gauge 
bosons of mass M. If the bosons are integrated out, then 
the resulting theory will be nonrenormalizable, i.e. it will 
be ill-behaved on momentum scales larger than the mass 
M. The nonrenormalizable theory corresponds to the 
order parameter LGW theory, except that in statistical 
mechanics one runs into infrared problems rather ultra- 
violet ones. Nevertheless, it turns out that in our case 
the critical behavior can still be determined exactly even 
after having integrated out the additional soft modes. 
The point is that the diffusive modes lead to an effec- 
tive long-range interaction between the order parameter 
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fluctuations that falls off in real space like 
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It is 



known that in general long-range interactions suppress 
fluctuation effects.H In our case they are strong enough 
to not only suppress quantum fluctuations, but also any 
remaining disorder fluctuations. The critical behavior is 
thus neither dominated by quantum fluctuations (since 
we work above the upper critical dimension dj ), nor by 
the disorder fluctuations, but rather is given by a simple, 
though non-standard (because of the long-range inter- 
actions) Gaussian theory. The resulting Gaussian fixed 
point allows for a correlation length exponent that satis- 
fies v > 2/d as required, and the exponents are dimen- 
sionality dependent for all d < 6. In d = 3 they are 
substantially different from either the mean-field expo- 
nents, or from those for a classical Heisenberg ferromag- 
net. This has striking observables consequences, as we 
will discuss. 

The outline of this paper is as follows. In Sec. |l| we 
first discuss some general aspects of itinerant ferromag- 
nets, and then we give our results for the critical expo- 
nents and for the equation of state near the critical point. 
Since the purpose of this paper is to give an exposition 
and discussion of these results that is as nontechnical as 
possible, they will be presented without any derivations. 



In Sec. Ill we discuss these results as well as several pos- 



sible experiments that could be performed to test our 
predicitions. Finally, in Sec. [j^ we sketch the derivation 
of our theoretical results. 



II. RESULTS 

In order to put the phase transition we are going to 
consider in perspective, let us first discuss the qualita- 
tive phase diagram that one expects for a disordered 
itinerant electron system in d = 3. Let F a < be the 




FIG. 1. Schematic phase diagram for a 3-d disordered itin- 
erant electron system in the plane spanned by the Landau 
parameter Fq and the disorder A at T — 0. See the text for 
further explanations. 



Fermi liquid parameter that characterizes the strength 
of the sytem's tendency towards ferromagnetism: For 
\Fq\ < 1 the system is paramagnetic with a spin sus- 
ceptibility x s - 1/(1 + F a ), while for \F$\ > I the clean 
Fermi liquid has a ferromagnetic ground state. In Fig. |l| 
we show the qualitative phase diagram one expects for a 
disordered system at T = in the Fff-X plane, where A is 
some dimensionless measure of the disorder. For A = 0, 
we have the transition from a paramagnetic metal (PM) 
to a ferromagnetic metal (FM) at Fq = — 1. At small but 
nonzero A this transition will occur at somewhat smaller 
values of \Fq |, since the disorder effectively increases the 
spin triplet electron-electron interaction amplitude, and 
hence \Fq\. This accounts for the downward curvature 
of the PM-FM transition line. At |F a | = 0, a metal- 
insulator transition of Anderson type is known to occur 
at a critical disorder value A c .E2l At nonzero \Fq | such a 
transition from a paramagnetic metal to a paramagnetic 
insulator (PI) still occurs, albeit it now is what is called 
an Anderson-Mott transition that takes place at a some- 
what larger value of the disorder. M The two transition 
lines will meet at a multicritical point M, and for large 
values of A and \Fq | one expects a ferromagnetic insu- 
lator (FI). The transitions from the FM and PI phases, 
respectively, to the FI phase have not been studied the- 
oretically, which is why we denote them by dashed lines 
in the figure. We will be mostly interested in the phase 
transition that occurs across the PM-FM transition line 
at finite disorder, but far away from the metal-insulator 
transition. However, in Sec. Ill below we will come back 



to the remaining regions in this phase diagram. 

In Fig. H we show the same phase diagram in the Fq - 
T plane for some value of the disorder < A << A c . 
With increasing temperature T, the critical value of \F§\ 
increases, since in order to achieve long-range order, a 
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FIG. 2. Schematic phase diagram for a disordered itinerant 
electron system in the plane spanned by the Landau parame- 
ter Fq and the temperature T. The inset shows the boundary 
of the critical region (dashed line) and the crossover line (dot- 
ted line) that separates classical critical behavior (cc) from 
quantum critical behavior (qc). 

larger \Fq\ is needed to compensate for the disorder- 
ing effect of the thermal fluctuations. The inset shows 
schematically the boundary of the critical region (dashed 
line) and the crossover line between classical and quan- 
tum critical behavior (dotted line). At any nonzero T, 
the asymptotic critical behavior is that of a classical 
Heisenberg magnet, but at sufficiently low T there is a 
sizeable region where quantum critical behavior can be 
observed. 

Our theoretical results for the zero temperature 
paramagnet-to-ferromagnet transition can be summa- 
rized as follows. Let t be the dimensionless distance from 
the line separating the regions PM and FM in Fig. [l]. 
Then the equation of state, which determines the mag- 
netization m as a function of t and the magnetic field h, 
can be written 



tm 



,d/2 . 



= h 



(2-1) 



where we have left out all prefactors of the various terms. 
Equation (2.1) is valid for all dimensions d > 2. Notice 
the term m d / 2 , which occurs in addition to what other- 
wise is an ordinary mean-field equation of state. It is a 
manifestation of the soft particle-hole excitations men- 
tioned in the Introduction. For d < 6 it dominates the 
m 3 -term, and hence we have for the exponent (3, which 
determines the vanishing of the zero-field magnetization 
via m(t,h= 0) ~ t?, 



f2/(d- 
P \ 1/2 



■2) for 2 < d < 6 
for d > 6 



(2.2a) 



Similarly, the exponent S, defined by m(t = 0, h) ~ h}/ s , 
is obtained as 
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for 2 < d < 6 
for d > 6 



(2.2b) 



Now let us consider the order parameter field Af (x, t) 
as a function of space and time, i.e. the field whose av- 
erage yields the magnetization, (M(x, t)) = m. Here 
the angular brackets (. . .) denote a trace with the full 
statistical operator, i.e. they include a quantum me- 
chanical expectation value, a disorder average, and at 
nonzero temperature also a thermal average. We first 
consider the case of T = 0, and Fourier transform to 
wave vectors q (with modulus q = |q|) and frequen- 
cies to. For the order parameter correlation function 
G(q,ui) = (M(q, u>) M(— q, — ui)) we find in the limit of 
small q and u), 



G(q,u) = 



1 



t + q d 2 + q 2 — icu/q 2 



(2.3) 



Here we have again omitted all prefactors of the terms in 
the denominator, since they are of no relevance for o ur 
discussion. The most interesting feature in Eq. ( |2.3| ) is 
the term q d ~ 2 . It is again an immediate consequence of 
the additi ona l soft mode s di scussed in the first section, 
and Eq. (|J), like Eq. (§t|), is valid for d > 2. For 
q = lu = 0, the correlation function G determines the 
magnetic susceptibility Xm ~ G(q = 0,u> = 0) in zero 
magnetic field. Hence we have Xm(t) ~ i" 1 ~ £~ 7 , where 
the last relation defines the critical exponent 7. This 
yields 



7=1 , 



(2.4) 



which is valid for all d > 2. 7 thus has its usual mean- 
field value. However, for nonzero q the anomalous q d ~ 2 
term dominates the usual q 2 dependence for all d < 4. 
The correlation function at zero frequency can then be 
written 



G(q,u = 0) 



1 



1 + («0 



d-2 



(2.5a) 



with the correlation length £ - t 1 /(<i-2) ^ t -v For d > 4 
the q 2 term is dominant, and we have for the correlation 
length exponent u, 



_/l/(d-2) 
1/2 



for 2 < d < 4 
for d > 4 



(2.5b) 



Note that v > 2 /d, as it must be according to the discus- 
sion in the Introduction. The wavenumber dependence 
of G at criticality, i.e. at t = 0, is characterize d by the 

2+)) 



exponent 77: 
obtain, 



G(q,cu = 0) ~ q~ 



From Eq. (2.3) we 



V = 



A-d for2<d<4 
. for d > 4 



(2.6) 



Finally, consider the correlation function at a wavenum- 
ber such that g£ = 1. Then it can be written 
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G(q = r 



1 



ILOT 



(2.7a) 



with the relaxation or correlation time r ~ £ 2 /t ~ 
where the last relatio n defi nes the dynami- 
cal critical exponent z. From Eq. ( 2.5t ) we thus obtain, 



d for 2 < d < 4 
4 for d > 4 



(2.7b) 



Notice that with increasing dimensionality d, the expo- 
nents v, r], and z 'lock into' their mean-field values at 
d = = 4, while /? and 5 do so only at d — d+ ++ = 6. 
In the special dimensions d = 4 and d = 6 the power law 
scaling behavior quoted above holds only up to additional 
multiplicative logarithmic dependences on the variables 
t, h, and T. Since these corrections to scaling occur only 
in unphysical dimensions they are of academic interest 
only, and we refer the interested reader to Refs. [| for 
details. 

The results for the clean case are qualitatively simi- 
lar, but the anomalous term in the equation of state, Eq. 
( |2.l[ ), is m d instead of m d / 2 . This is because the ad- 
ditional soft modes in that case are ballistic instead of 
diffusive, so their frequency scales with wavenumber like 
to ~ q rather than to ~ q 2 . As a result, the two special di- 
mensions d+ + and coincide, and are now d+ + = 3, 
while the upper critical dimension proper, above which 
fluctuations are irrelevant, is g?+ = 1. For 1 < d < 3, 
the exponent values are = u = l/(d— 2), 6 = z = d, 
rj = 3 — d, and 7 = 1. For d > 3, all exponents take on 
their mean-field values as they do in the disordered case 
for d > 6, and in d = 3 there are logarithmic corrections 
to power-law scaling. 

We now turn to the behavior at nonzero temperatures. 
Then the equation of state acquires temperature correc- 
tions, and it is helpful to distinguish between the cases 
m » T and m << T, with m and T measured in suit- 
able units. Taking into account the leading corrections 
in either limit, the equation of state reads 

tm + m d/2 (1 + T/m) = h (for m»T) , 



(t + T( d - 2 )/ 2 ) m + ■ 



(for T»m) . (2. 



Equation (2.8) shows that for any nonzero temperature 



the asymptotic critical behavior is not given by the quan- 
tum critical exponents. Since Eq. fl2.8| ) takes temper- 
ature into account only perturbatively, it correctly de- 
scribes only the initial deviation from the quantum crit- 
ical behavior, and approximates the classical critical be- 
havior by the mean-field result. A full crossover calcula- 
tion would yield instead the classical Heisenberg critical 
behavior in the asymptotic limit. Also, we are consid- 
ering only the saddle point contribution to the magneti- 
zation. For models with no additional soft modes it has 
been shown that fluctuations that act as dangerous irrel- 
evant variables introduce another temperature scale-that 
dominates the one obtained from the saddle point .BEJ In 
the present case, however, fluctuations are suppressed by 



the long-range nature of the effective field theory, and 
the fluctuation temperature scale is subdominant. The 
behavior described by Eq. (2.8) can be summarized by 



means of a generalized homogeneity law, 

m(t, T, H) = b- p f v m(tb x l\ Tb*' 1 ' ', Hb Sfi/u ) . (2.9a) 

Here 0, v, and 5 have the values given above, and b is an 
arbitrary scale factor. 



2v 



(2.9b) 



is the crossover exponent that describes the deviation 
from the quantum critical behavior due to the relevant 
perturbation provided by the nonzero temperature. The 



entry Tb^ v = Tb 2 in the scaling function in Eq. fl2.9a| ) 
reflects the fact that the temperature dependence of the 
saddle point solution is determined by that of the dif- 
fusive modes, i.e. frequency or temperature scales like 
T ~ q 2 ~ b~ 2 . The critical temperature scale, T ~ b~ z , 
would be dominant if it were present, but since the lead- 
ing behavior of the magnetization is not determined by 
critical fluctuations, it is suppressed. 

By differentiating Eq. (2.9a) with respect to the mag- 
netic field h, one obtains an analogous homogeneity law 
for the magnetic susceptibility, \ m , 

Xm (t, T, H) = tfl v Xmitb 1 ^, Tb*/», Hb & PI») , (2.10a) 

with 



7 = /3(<5-l) = l , 



(2.10b) 



in agreement with Eq. (2^). This result is in agreement 
with a more direct calculation of x-m- The same tem- 
perature corrections that modify the equation of state, 
Eq. (2.8), lead to a replacement of the term q d ~ 2 in the 
denominator of Eq. ( O ) by (q 2 + T)^ 2 )/ 2 . Since the 
homogeneous order parameter correlation function de- 
termines the spin or order parameter susceptibility, this 
yields 



Xm(t, T) 



1 



(2.10c) 



in agreement with Eqs. (2.10a, 2.10h| ). 

Finally, the critical behavior of the specific heat cv 
has been calculated. It is most convenient to discuss the 
specific heat coefficient, = lim-j^g c v /T, which in a 
Fermi liquid would simply be a constant. Its behavior at 
criticality, t — 0, is adequately represented by the integral 



dq 



T + q d + qi + /l 1 - 1 ^ 



(2.11a) 



Remarkably, in zero magnetic field, jv diverges logarith- 
mically as T — > for all dimensions 2 < d < 4. This can 
be shown to be a consequence of the dynamical expo- 
nent z being exactly equal to the spatial dimensionality 
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d in that range of dimensionalities. If one restores the 
dependence of 7y on t, then one obtains a generalized 
homogeneity law with a logarithmic correction for the 
leading scaling behavior of 7y, 



j v (t,T,H) = 6(4- d) hxb 

+F 1 {tb 1/v ,Tb z ,Hb S0/u ) 



(2.11b) 



Here Q(x) denotes the step function, and F 7 is an un- 
known scaling function. Note that jv is determined by 
Gaussian fluctuations and depends on th e crit ical tem- 
perature scale, i.e. T scales like t vz in Eq. ( [2.11b ). This is 
the leading temperature scale, and whenever it is present 
it dominates the diffus ive temperature scale that shows 
in Eqs. (|J) and ( p-10D - 

In the clean case, Eqs. ( 2.9a ) and ( [2.10 ) still hold, if 
one uses t he appropriate expon ent va lues and replaces 
Eq. j2.9bD by <j) = v. In Eq. ( |2.11a| ), the term q 4 in 
the denominator of the integrand gets replaced by q 3 , 
and consequently the argument of the 9-function in Eq. 
( p.llbj ) is 3 - d rather than 4 - d. 



III. EXPERIMENTAL IMPLICATIONS, AND 
DISCUSSION 

A. Experimental Implications 

Let us now discuss the experimental implications of the 
results presented in the preceding section. Obviously, one 
needs a material that shows a transition from a paramag- 
netic state to a ferromagnetic one at zero temperature as 
a function of some experimentally tunable parameter x. 
Obvious candidates are magnetic alloys of the stoichiom- 
etry P^Fi-^, with P a paramagnetic metal and F a ferro- 
magnetic one. Such materials show the desired transition 
as a function of the composition parameter x; examples 
include Ni for the ferromagnetic component, and Al or 
Ga for the paramagnetic one.E3 At the critical concen- 
tration x c they also are substantially disordered, but due 
to the fact that both constituents are metals they are far 
from any metal-insulator transition. Our theory should 
therefore be applicable to these systems. The schematic 
phase diagram at T = in the T-x plane is shown in Fig. 
||. Notice that this is a realistic phase diagram, as op- 
posed to the 'theoretical' ones in Figs. [l]and|^. A change 
of the composition parameter x leads, besides a change of 
Fq , to many other changes in the microscopic parameters 
of the system. As x is varied, the system will therefore 
move on a complicated path in the diagram shown in, 
say, Fig. [|. However, since the critical behavior near the 
transition is universal, it is independent of the exact path 
traveled. 

One possible experiment would consist in driving the 
system at a low, fixed temperature through the transi- 
tion by changing the composition x. While this involves 
the preparation of many samples, this way of probing 




FIG. 3. Schematic phase diagram for an alloy of the form 
Pi Fi- X . T c is the Curie temperature for the pure ferromag- 
net F, and x c is the critical concentration. 



a quantum phase transition has been used Jjo observe 
the metal-insulator transition in P doped Si. til It might 
also be possible to use the stress tuning technique that 
has been used for the same purpose.Ej Either way one 
will cross the transition line along a more or less vertical 
path in Fig. |[ and for a sufficiently low temperature this 
path will go through both the classical and the quantum 
critical region indicated in the inset in Fig. |^. Due to 
the large difference between the quantum critical expo- 
nents quoted in Sec. [n| and the corresponding exponents 
for classical Heisenberg magnets, the resulting crossover 
should be very pronounced and e asily observable. For 
instance, for 3-d systems our Eq. ( 2.2a ) predicts (3 — 2, 
while the value for the thermal transition is (3 c i ass « 0.37. 
The resulting crossover in the critical behavior of the 
magnetization is schematically shown in Fig. || Alterna- 
tively, one could prepare a sample with a value of x that 
is as close as possible to x c , and measure the magnetic 
field dependence of the magnetization, extrapolated to 
T = 0, to obtain the exponent 8. Again, there is a large 
difference between our prediction of S — 1.5 in d = 3, 
and the classical value <5 c iass ~ 4.86. 

Yet another possibility is to measure the zero-field 
magnetic susceptibility as a function of both t — \x — x c \ 
and T. Equation ( |2.10a ) predicts 



Xm (t,T)=T 



(3.1) 



Here f x is a scaling function that has two branches, 
for x > x c , and /~ for x < x c . Both branches ap- 
proach a constant for large values of their argument, 
fx[y — > oo ) = const. For small arguments, we have 
fx(y ~~ * 0) ~ \/y, while /~ diverges at a nonzero value 
y* of its argument that signalizes the classical transition, 

fx(y y*) ~ (y - y*r 7clas % with 7class « 1.39 the 

susceptibility exponent for the classical transition. Our 
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FIG. 4. Schematic critical behavior of the magnetization 
m at nonzero temperature, showing the crossover from the 
quantum critical behavior (/3 = 2, dashed line) to the classical 
critical behavior (/? ps 0.37, dotted line). Notice that the 
actual transition is classical in nature. 

prediction is then that a plot of \m T 1 ^ 2 versus T/t 2 will 
yield a universal function the shape of which is schemat- 
ically shown in Fig. |[ Notice that the exponents are 
known exactly, so the only adjustable parameter for plot- 
ting experimental data will be the position of the critical 
point. This is on sharp contrast to some other quantum 
phase transitions, especially metal- insulator transitions, 
where the exponent values are not even approximately 
known, which makes scaling plots almost meaningless. Ej 
Finally, one can consider the low-tempe rature behav- 
ior of the specific heat. According to Eq. ( |2.11b| ), as the 
temperature is lowered for x > x c the leading tempera- 
ture dependence of the specific heat will be 

c v {T)~T\nT . (3.2a) 

At criticality this behavior will continue to T = 0, while 
for x > x c it will cross over to 

c y (T) - (kit) T . (3.2b) 

For x < x c one will encounter the classical Heisenberg 
transition where the specific heat shows a finite cusp (i.e., 
the exponent a, defined by cy ~ (T— T c )~", is negative). 

B. Theoretical Discussion 

There are also various theoretical implications of the 
results presented in Sec. Q. One aspect is the general 
message that the usual LGW philosophy must not be ap- 
plied uncritically to quantum phase transitions, because 
of the large number of soft modes that exist at zero tem- 
perature in a generic system. If any of these couple to 
the order parameter, then an effective theory entirely in 
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FIG. 5. Schematic prediction for a scaling plot of the mag- 
netic suscteptibility. 

terms of the order parameter will not be well behaved. In 
the present case we have actually been able to use this to 
our advantage, since the long-ranged interaction that the 
additional soft modes induce in the order parameter the- 
ory suppress the disorder fluctuations, which is the reason 
for the remarkable fact that we are able to exactly deter- 
mine the critical behavior of a three-dimensional, disor- 
dered system. In general, however, the presence of soft 
modes in addition to the order parameter fluctuations 
will call for the derivation of a more complete low-energy 
effective theory that keeps all of the soft modes explicitly. 

Another very interesting aspect is a connection be- 
tween our results on the ferromagnetic transition, and 
a substantial body of literature on a problem that ap- 
pears in the theory of the metal-insulator transition in 
interacting disordered electron systems, i.e. the transi- 
tion from PM to PI in Fig. [p. This problem has been 
known ever since the metal-insulator transition of inter- 
acting disordered electrons was first considered, and it 
has led to substantial confusion in that field. Early work 
on the metal-insulator transition showed that in two- 
dimensional systems without impurity spin-flip scatter- 
ing, the spin-triplet interaction amplitude scaletLto large 
values under renormalization group iterations O This is 
still true in d = 2 + e, and since the runaway flow oc- 
curs before the metal-insulator transition is reached, this 
precluded the theoretical description of the latter in such 
systems. This problem was interpreted, incorrectly as it 
turned out later, as a signature of local moment forma- 
tion in all dimensions .113 Subsequently, the j-aresent au- 
thors studied this problem in some detailB We were 
able to explicitly resum the perturbation theory and show 
that at a critical value of the interaction strength, or of 
the disorder, there is a bulk, thermodynamic phase tran- 
sition in d > 2 that is not the metal-insulator transition. 
While this ruled out local moments (which would not 



7 



lead to to phase transition) , the physical meaning of this 
transition was obscure at the time since no order param- 
eter had been identified, and its description was entirely 
in terms of soft diffusion modes. However, the critical 
exponents obtained are identical to those given in Sec. H 
for the quantum ferromagnetic phase transition, and in, 
both cases logarithmic corrections to scaling are found. E2l 
Because the exponents in the two cases are identical, we 
conclude that the transition found earlier by us, whose 
physical nature was unclear, is actually the ferromagnetic 
transition. One also concludes that our speculations in 
Rcf. [11] about the nature of the ordered phase as an 'in- 
completely frozen spin phase' with no long-range mag- 
netic order, were not correct; this phase is actually the 
metallic ferromagnetic phase. On the other hand, the 
techniques used in that reference allowed for a determi- 
nation of the qualitative phase diagram as a function of 
dimensionality which our present analysis is not capable 
of. This analysis showed the existence of yet another in- 
teresting dimensionality above d = 2, which we denote 
by d* . With the appropriate reinterpretation of the 'in- 
completely frozen spin phase' as the ferromagnetic phase, 
the qualitative phase diagram for 2 < d < d* is shown in 
Fig. ^]. Compared to Fig. [l], the following happens as d 
is lowered from d = 3: The multicritical point M moves 
to downward, and at d = d* it reaches the A-axis. d* was 
estimated in Ref. [l^ to be approximately d* = 2.03. For 
d < d*, the insulator phase can therefore not be reached 
directly from the paramagnetic metal. This explains why 
in the perturbative renormalization group calculations 
in d = 2 + e one necessarily encounters the ferromag- 
netic transition first, and it should finally put to rest the 
long discussion about the physical meaning of the run- 
away flow that is encountered in these theories. It also 
shows that none of these theories are suitable for studying 
the metal-insulator transition in the absence of spin-flip 
mechanisms, as they start out in the wrong phase. 

It should also be pointed out that our earlier theory 
depended crucially on there being electronic spin con- 
servation. This feature would be lost of there were some 
type of impurity spin flip scattering process. In that case, 
the soft modes that lead to the long-range order param- 
eter interactions acquire a mass or energy gap, and at 
sufficiently large scales the interactions are effectively of 
short range. The asymptotic critical phenomena in this 
case are described by a short-range, local order parameter 
field theory with a random mass, or temperature, term. 
Such a term is present in the case of a conserved order 
parameter also, but due to the long ranged interaction it 
turns out to be irrelevant with respect to the nontrivial 
Gaussian fixed point. In the absence of the conservation 
law, however, the random mass term is relevant with re- 
spect to the Gaussian fixed point analogous to the one 
discussed here. This underscores the important role that 
is played by the order parameter being conserved in our 
model. The quantum phase transition in a model where 
it is not has been discussed in Ref. ||. 

We finally discuss why some of our results are in dis- 
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FIG. 6. Schematic phase diagram for a disordered itiner- 
ant electron system at T = close to d = 2. The phases 
shown are the paramagnetic metal (PM), the ferromagnetic 
metal (FM), and the insulator (I) phase. Whether within I 
there is another phase transition from a ferromagnetic to a 
paramagnetic insulator is not known. 

agreement with Sachdev'si general scaling analysis of 
quantum phase transitions with conserved o rder param - 
eters. For instance, it follows from our Eqs. ( pC[ , [Dig ) 
that the Wilson ratio, defined as W = (m/H)/(C v /T), 
diverges at criticality rather than being a universal num- 
ber as predicted in Ref. |[ Also, for 2 < d < 4 the 
function F 7 in Eq. ( |2.11b| ), for t — and neglecting cor- 
rections to scaling, is a function of T / in agreement 
with Ref. |6| but for d > 4 this is not the case. The 
reason for this breakdown of general scaling is that we 
work above an upper critical dimensionality, and hence 
dangerous irrelevant variables^ appear that prevent a 
straightforward application of the results of Ref. [sj to the 
present problem. These dangerous irrelevant variables 
have to be considered very carefully, and on a case by 
case basis. This caveat is particularly relevant for quan- 
tum phase transitions since they tend to have a low upper 
critical dimension. It is well known that a given irrelevant 
variable can be dangerous with respect to some observ- 
ables but not with respect to others. Specifically in our 
case there is a dangerous irrelevant variable that affects 
the leading scaling behavior of the magnetization, but 
not that of the specific heat coefficient, and this leads to 
the divergence of the Wilson ratio. This dangerous ir- 
relevant variable is also the reason why the exponents (3 
and <5, which describe the critical behavior of the magne- 
tization, remain dimensionality dependent up to d — 6, 
while all other exponents 'lock into' their mean-field val- 
ues already at d = 4. 
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IV. THEORETICAL OUTLINE 

Here we sketch the derivation of the results that were 
presented in Sec. |l| We do so for completeness only, and 
will be very brief. A detailed account of the derivation 
can be found in Ref. ||for the disordered case, and in Ref. 
H for the. clean case. 

Hcrtza has shown how to derive an LGW functional 
for a quantum ferromagnet. One starts by separating the 
spin-triplet part of the electron-electron interaction, i.e. 
the interaction between spin density fluctuations, from 
the rest of the action, writing 



S — S + SjL 



with 



S int = y / dxn s (x) -n s (a 



(4.1a) 



(4.1b) 



Here S^ t is the spin-triplet interaction part of the ac- 
tion, and 5*0 contains all other parts, in particular the 
electron-electron interaction in all other channels. r t 
is the spin triplet interaction amplitude, which is re- 
lated to the Landau paramter Fq used above by T t — 
—Fq/(1 + Fq), n s (x) is the electron spin density vec- 
tor, x = (x, r) denotes space and imaginary time, and 

J dx — f dx Jq^ T dr. In the critical region near a quan- 
tum phase transition, imaginary time scale like a length 
to the power z, and the space-time nature of the integrals 
in the action accounts for the system's effective dimen- 
sion d + z that was mentioned in the Introduction. 



Now S, 



J int i s decoupled 
Stratonovich transformat 



ion 1 ! 



means of a Hubbard- 
The partition function, 
apart from a noncritical multiplicative constant, can then 
be written 



D[M] exp(-$[M]) 



(4.2a) 



with the LGW functional 

r, 



In yexp 



dx M(x) ■ M(x) 
-T t J dx M(x) ■ n s (x) 



(4.2b) 



So 



Here (. . .)<?„ denotes an average taken with the action Sq. 
If the LGW functional is formally expanded in powers 
of M, then the term of order M™ obviously has a coef- 
ficient that is given by a connected n-point spin density 
correlation function of the 'reference system' defined by 
the action So- 

At this point we need to remember that our reference 
system So contains quenched disorder, which has not 
been averaged over yet. The n-point correlation func- 
tions that form the coefficients of the LGW functional 
therefore still depend explicitly on the particular realiza- 
tion of the randomness in the system. The average over 



the quenched disorder, which we denote by {. . .}dis, re- 
quires averaging the free energy, i.e. we are interested 
in {In Z}dis- r-This is most easily done by means of the 
replica trick,E3, i.e. one writes 



lim - [{Z"}dis - 1] 

n— >0 n 



lim — 

n— >0 n 



{lnZ} d 



- 1 



(4.3) 



where the index a labels n identical replicas of the sys- 
tem. The disorder average is now easily performed by 
expanding the exponential in Eq. (4.2). Upon reexpo- 
nentiation, the coefficients in the replicated LGW func- 
tional are disorder averaged correlation functions of the 
reference system that are cumulants with respect to the 
disorder average. The Gaussian part of $" is simply 



$f 2) [M Q ] = IJ d Xl dx 2 M Q (xi) 



S(xi 



X2) 



-T t x(xi - x 2 ) 



■M a (x 2 ) 



(4.4) 



Here x(x) is the disorder averaged spin susceptibility or 
2-point spin density correlation function of the reference 
system. The cubic term, <I>?U has a coefficient given by 
the averaged 3-point spin density correlation function. 
For the quartic term, the cumulant nature of these corre- 
lation functions leads to two terms with different replica 
structures, and higher order terms have correspondingly 
more complicated structures. 

The next step is to calculate the spin density corre- 
lation functions for the reference system. It now be- 
comes important that we have kept in our action Sq 
the electron-electron interaction in all channels except for 
the spin-triplet one that has been decoupled in deriving 
the LGW functional. At this point our treatment devi- 
ates from that of Hertz, who took the reference ensemble 
to describe noninteracting electrons. This was generally 
considered an innocent approximation that should not 
have any qualitative effects. However, this belief was mis- 
taken, since the spin density correlations of interacting 
electrons are qualitatively different from those of nonin- 
teracting ones. The spin susceptibility can be easily cal- 
culated in perturbation theory. The result shows that the 
static spin susceptibility as a function of the wavenumbcr 
q is nonanalytic at q = 0. For small q it has the form 



x(q) = const 



(4.5) 



The nonanalyticity is a consequence of the presence of 
soft particle-hole excitations in the spin-triplet channel, 
and it occurs only in an interacting electron system. That 
is, the prefactor of the q d ~ 2 term, which we have sup- 



pressed in Eq. (4.5), vanishes for vanishing interaction 
amplitudes. Rcnormalization group arguments can then 
be used to ascertain that this perturbative result indeed 
represents the exact behavior of % in the long- wavelength 
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limit. If one also considers the frequency dependence of 
X, one obtains the Gaussian part of the LGW functional 
in the form 



$ f2)[M]=^^M Q (q, W „) [t Q + q 



d-2 



where 



-q 2 + \uj n \/q 2 ] ■M a (-q,,-LJ n ) 



t Q = l-T t Xs (q^0,cj„ = 0) 



(4.6a) 



(4.6b) 



is the bare distance from the critical point, and the uo n = 
2nTn are bosonic Matsubara frequ encies. 

The Gaussian theory, Eqs. ( |4.6| ), can be analyzed us- 
ing standard renormalization group techniques B3 Such 
an analysis reveals the existence of a Gaussian fixed point 
whose critical properties are the ones given in Sec. ||. 
The remanining question is whether this fixed point is 
stable with respect to the higher, non-Gaussian terms in 
the action. These higher terms also need to be considered 
in order to obtain the critical behavior of the magnetiza- 
tion. 

A calculation of the higher correlation functions that 
determine the non-Gaussian vertices of the field theory 
shows that the nonanalyticity that is analogous to the 
one in the spin susceptibility, Eq. (4.5), is stronger and 
results in a divergence of these correlation functions in 
the zero frequency, long- wavelength limit. Specifically, 
the leading behavior of the n-point spin density correla- 
tion that determines the coefficient of the term of order 
M™ in the LGW functional, considered at vanishing fre- 
quencies as a function of a representative wavenumber q, 



X (n \q^0) 



d+2-2n 



(4.7) 



As a result, the coefficients cannot, as usual, be expanded 
about zero wavenumber, and the theory is nonlocal. De- 
spite this unpleasant behavior of the field theory, it is 
easy to see by power counting that all of these terms ex- 
cept for one are irrelevant with respect to the Gaussian 
fixed point in all dimensions d > 2. The one exception 
is the quartic cumulant contribution that is the disorder 
average of the square of the spin susceptibility, which is 
marginal in d = 4, but irrelevant in all other dimensions. 
This term is physically of interest, since it represents the 
random mass or random temperature contribution that 
one would expect in a theory of disordered magnets, and 
that was mentioned in Sec. Ill B| above. 

The conclusion from these considerations is that apart 
from logarithmic corrections to scaling in certain special 
dimensions, the Gaussian theory yields the exact criti- 
cal behavior, and the only remaining question pertains 
to the form of the equation of state. Since the quar- 
tic coefficient x^ is a dangerous irrelevant variable for 
the magnetization, this requires a scaling interpretation 
of the infrared divergence of x (4) ■ In Ref. | it has been 



show n that for scaling purposes the wavenumber q in Eq. 
(4.7) can be identified with the magnetization m 1 / 2 . This 
is physically plausible, since the divergence stems from an 
integration over soft modes that are rendered massive by 
an external magnetic field. Since a nonzero magnetiza- 
tion acts physicall y lik e a magnetic field, it cuts off the 
singularity in Eq. (4.7). With this interpretation of the 
singular coefficients, the term of order m n in the sad- 
dle point solution of the LGW theory has the structure 

which leads to the equation 



l/2\d+2-2n 



T d/2 



of state as given in Eq. (2.1). One might wonder why the 
magnetic fluctuations in the paramagnetic phase do not 
also cut off the singularity in Eq. (4.7), and thus weaken 
or even destroy the effects discussed above. While such 
a cutoff mechanism does exist, it enters the theory only 
via the fluctuations, which are RG irrelevant with respect 
to the Gaussian fixed point. It therefore shows only in 
the corrections to scaling, but not in the leading critical 
behavior. 

Again, all of these arguments can be repeated for the 
case without disorder. The only changes one encounters 
pertain to the values of various exponents due to the 
different character of the soft modes. This leads to the 
results quoted in Sec. Q. 
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